Absolute values of the London penetration depth in YBa2Cu306+y measured by zero 
field ESR spectroscopy on Gd doped single crystals 
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Zero-field electron spin resonance (ESR) of dilute Gd ions substituted for Y in the cuprate super- 
conductor YBa2Cu306-i-y is used as a novel technique for measuring the absolute value of the low 
temperature magnetic penetration depth A(r — > 0). The Gd ESR spectrum of samples with « 1% 
substitution was obtained with a broadband microwave technique that measures power absorption 
bolometrically from 0.5 GHz to 21 GHz. This ESR spectrum is determined by the crystal field that 
lifts the level degeneracy of the spin 7/2 Gd''"'" ion and details of this spectrum provide information 
concerning oxygen ordering in the samples. The magnetic penetration depth is obtained by relating 
the number of Gd ions exposed to the microwave magnetic field to the frequency-integrated inten- 
sity of the observed ESR transitions. This technique has allowed us to determine precise values 
of A for screening currents flowing in the three crystallographic orientations (a, h and c) in sam- 
ples of GdxYi_xBa2Cu306+y of three different oxygen contents y = 0.993 (Tc = 89 K), y = 0.77 
{Tc = 75 K) and y = 0.52 (Tc = 56 K). The in-plane values are found to depart substantially from 
the widely reported relation Tc oc 1/ . 



I. INTRODUCTION 

The magnetic penetration depth A figures prominently 
in many aspects of superconductivity. It is the length 
scale over which an external magnetic field is screened by 
a superconductor in the Meissner state, it sets the size of 
vortices in the mixed state, and it is a controlling factor in 
many other properties such as microwave and far infrared 
absorption. The fundamental importance of the penetra- 
tion depth is that it provides a rather direct measure of 
the superfluid phase stiffness ps oc or what is of- 

ten referred to as superfluid density ris/m* = l//xoe^A^, 
which can be a tensor quantity due to effects such as an 
anisotropic effective mass m*. In the cuprates, the de- 
pendence of A on temperature T and hole-doping have 
provided major pieces of the cuprate puzzle. Microwave 
measurements of AA(T) at low T revealed a linear de- 
pendence providing some of the first key evidence of the 
presence of line nodes in the superconducting gap in the 
cuprate^' and the temperature dependence of A(T) near 
Tc suggested that the transition is governed by 3DXY 
critical fluctuationsi^i Evidence that the low temperature 
value of the phase stiffness Ps{T — > 0) depends linearly 
on Tc, has been obtained from muon spin relaxation mea- 
surements of A(T)i^ 

Accurate and precise measurements of the absolute 
value of A(T — > 0) are particularly important for un- 
derstanding superconductivity in the cuprates, and are 
critical for broad studies of doping dependence. Having 
an absolute value of A is also crucial in the interpretation 
of other experiments that involve superfluid screening, 
with examples including the microwave conductivity, the 
lower critical field Hci, and the conversion of measure- 
ments of AA(T) to the temperature dependent superfluid 
density cx 1/A^(T). However, the absolute value of A has 
proven to be particularly difficult to measure. A wide va- 



riety of techniques exist for doing so, but each has impor- 
tant drawbacks that must be acknowledged. One class 
of techniques takes place in the mixed state, probing the 
length scale over which fields decay away from a magnetic 
vortexi4 Muon spin relaxation {p,SR) falls into this class 
and suffers from the difficulty that the screening length 
scale obtained from a mixed state measurement in rela- 
tively high fields can differ from the penetration depth 
in the Meissner state because of non-linear and non-local 
effects.^ Vortex state measurements are also unable to di- 
rectly resolve the anisotropy in A in orthorhombic materi- 
als because the screening currents circulate around a vor- 
tex. A second class of techniques measures magnetic fiux 
exclusion from a sample in the Meissner state, essentially 
by measuring the magnetic susceptibility at DC, radio or 
microwave frequencies These can precisely determine 
the temperature dependence AA(T) = A(T) — A(T'o) rela- 
tive to a base temperature Tg. However, a measurement 
of the absolute value amounts to a comparison between 
the physical volume of the sample and a susceptibility de- 
termination of the volume that is field-free, the difference 
being the small volume at the surface that is penetrated 
by the field. Such a measurement is most sensitive if the 
ratio of surface area to volume is large, which can be 
achieved by working with powders. However, difficulties 
in aligning powders and modelling the shape and size dis- 
tribution of grains introduce substantial uncertainties in 
the absolute value of Xfi For a macroscopic, high quality 
crystal of a cuprate superconductor, it is nearly impos- 
sible to measure sample dimensions precisely enough to 
achieve an absolute measurement of A this way, except 
in geometries dominated by the very large values of A for 
screening currents running in the c-direction. 

A very direct measurement of the penetration depth 
has been achieved by measuring the transmission of low 
frequency fields through very thin films using mutual 
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inductance techniques. While this is a very sensitive 
method, its chief drawbacks are the need to work with 
films and the difficulty of measuring anisotropyiS Tech- 
niques that can resolve the anisotropy include far infrared 
and optical measurements that determine the inductive 
response of the superfluid at high frequencies piS*ii Such 
measurements are difhcult because of their need for spec- 
tra over a very wide frequency range in order to per- 
form Kramers-Kronig transforms and are also suscepti- 
ble to error because the measured response can also in- 
clude screening by uncondensed carriers in addition to 
the superfluid. However, this approach can be used to 
determine the anisotropy of A simply by polarizing the 
light and measuring optical properties along all of the 
principal axes. Such infrared determinations of the ab- 
solute value of A at low temperature have been com- 
bined with microwave measurements of AA(r) to pro- 
duce the anisotropic 1/A^(T) in YBa2Cu306+y at y=0.60 
and 0.95,^^ but it would be far preferable to have a more 
direct microwave measurement that gives the absolute 
value of A, measured on the same sample as that used to 
determine the temperature dependence. 

This paper describes a microwave technique that di- 
rectly measures the small volume near the surface that 
is field-penetrated in the Meissner state, rather than in- 
ferring A from a measurement of the very large volume 
that is field-free (see Fig. P). We use this technique here 
to measure the absolute value of A for all three princi- 
pal axes in YBa2Cu306+y. The key to the method is to 
embed a randomly distributed low concentration of mag- 
netic moments in the superconductor in such a way that 
they will act as local field probes and yet will minimally 
affect the transport properties. The magnetic moments 
will serve as a non-interacting spin system whose electron 
spin resonance (ESR) spectrum is a sum of one-particle 
ESR transitions. The energy levels of a single spin are 
determined by its crystallographic environment, i.e. the 
charge distribution around the magnetic ion, and the re- 
sulting susceptibility is that of a single spin multiplied by 
the number of spins that participate in the ESR process. 
If the moments are evenly distributed in the sample, the 
number of spins that participate in the process is pro- 
portional to the fraction of the sample's volume that is 
exposed to the field, i.e. the volume within a penetra- 
tion depth of the surface. In order to count the number 
of spins exposed, we measure the zero field ESR absorp- 
tion spectrum with a microwave frequency magnetic field 
applied perpendicular to the principal spin axis acting as 
a small perturbation. 

In the next section of this paper we describe the sam- 
ples, the broadband microwave apparatus and the mea- 
surements of the apparent microwave surface resistance 
that ultimately yield values of the penetration depth. 
The third section describes how the theoretical ESR spec- 
trum was calculated using an effective crystal field Hamil- 
tonian, which is a necessary prerequisite for extracting A. 
In the fourth section we examine the issue of oxygen or- 
dering in the CuO chains of YBa2Cu306-(-y and comment 




FIG. 1: Schematic drawing of a platelet sample of a supercon- 
ductor having an anisotropic penetration depth A in a uniform 
applied microwave magnetic field Hrf. In the above picture, 
the field penetrates a volume determined by Xa and Ac and 
the crystal dimensions £a and Ic, but contains no contribution 
from Afc. By combining measurements where we rotate and 
cleave the crystal, we are able to determine values of A for all 
three crystallographic directions. 

upon the sensitivity of our measurements to this order- 
ing. Finally, we summarize our results and present new 
values for A in all three directions, for three different dop- 
ings in the YBa2Cu306+y system. These results indicate 
a substantial departure from Tc oc 1/A^, the so-called 
Uemura scaling. 

II. SAMPLE PREPARATION AND 
EXPERIMENTAL TECHNIQUES 

The objective of this work is to provide accurate mea- 
surements of the intrinsic absolute value of the mag- 
netic penetration depth in the limit T —^ 0. It is well- 
known that impurity doping can have a strong influ- 
ence on the low energy density of states in a d-wave 
superconductor, and for the present work care must 
be taken to ensure that the introduction of the spin-probe 
impurity does not alter A(r —> 0). This constrains the 
Gd concentration x, but conversely we wish to maximize 
the amplitude of the ESR response, which scales as x. 
We have found a good compromise at a nominal value of 
x« 1% which results in an easily resolved ESR spectrum 
and does not significantly alter the intrinsic low temper- 
ature properties of the material. 

In GdxYi_xBa2Cu306-i-y , the dopant Gd'^^ ion sub- 
stitutes for the Y^+ ion which is sandwiched between 
the two Cu02 planes in the YBa2Cu306+y unit cell that 
support the bulk of the electronic transport. The energy 
cost for cross-substitution of the Gd on the Ba site is pro- 
hibitive and such cation substitution occurs only for the 
lighter and larger members of the rare earth series. The 
ESR spectrum is determined completely by the crystal 
field environment of the Gd'^^ ion and therefore know- 
ing that there is only one site for substitution means that 
changes in the ESR spectrum from sample-to-sample can 
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only be due to changes of the oxygen order in the CuO 
chains. Another important feature of substitution on the 
Y site is that the Gd impurities are expected to be weak 
(Born) hmit quasiparticle scattering centers as was found 
to be the case for Ca^+ substitution on the same site»i^ 
Born scatterers have httle effect on the superfluid den- 
sity in the low temperature limit. This is in contrast, for 
instance, to the case of the non-magnetic impurity TjV?^ 
which substitutes into the Cu02 plane, and is known to 
be strongly pair-breaking. 

The GdxYi_xBa2Cu306-Hy samples were grown us- 
ing a self-flux method in BaZrOs crucibles as de- 
scribed elsewhereiifii The starting materials had at least 
99.999 at. % purity with the nominal concentration of 
Gd ions set by the ratio of oxide precursor substitution 
Gd2 03:Y2 03. The as-grown platelet single crystals were 
mechanically detwinned, annealed in flowing oxygen to 
set the oxygen content, and given a final anneal at lower 
temperature for oxygen ordering. The annealing parame- 
ters have been previously published for the slightly over- 
doped (nearly full CuO chains, y=0.993)^''^ and ortho- 
II ordered (alternating full and empty chains, y=0.52)^^ 
crystals. Other ordered phases at intermediate doping 
can also be produced with larger chain-ordered super- 
lattices, namely ortho-III^i and ortho-VIII.^^ Crystals 
at these doping levels were produced using procedures 
very similar to that used to produce the ortho-II samples, 
however the initial temperatures for setting the oxygen 
content were different (668 °C for ortho-VIII and 609 °C 
for ortho-III) as were the final low temperature anneal- 
ing temperatures used to establish oxygen ordering in the 
CuO chain layers (35 °C for ortho-VIII (y=0.67) order 
and 60 °C for ortho-III (y=0.77) order). 

A critical parameter in our data analysis is the ab- 
solute concentration of Gd ions x, which was measured 
using electron-probe micro-analysis (EPMA). We used a 
fully- automated CAMECA SX-50 instrument, operated 
in the wavelength-dispersion mode with the following op- 
erating conditions: excitation voltage, 15 kV; beam cur- 
rent, 20 nA; peak count time (240 s for Gd), 20 s; back- 
ground count-time, 10 s each side of peak (120 s for Gd); 
spot diameter, 5 /xm. Data reduction was done using the 
"PAP" 0(/9Z) method.^" For the elements considered, the 
following standards, X-ray lines and crystals were used: 
GagGdgOia, GdLa, LiF; YBa2Cu3 06.95, BaLa, PET; 
YBa2Cu3 6.95, CuKa, LiF; YBa2Cu3 06.95, OKa, W/Si 
multilayer dispersion element. For the present study, we 
found it to be important that each GdxYi_xBa2Cu306+y 
crystal was examined individually. Because the sample 
must be embedded in epoxy for EPMA, this was done fol- 
lowing all microwave measurements. For a nominal Gd 
concentration of x= 1% in the growth flux, the resulting 
concentration in the crystals was found vary to over the 
range from x= 1.05(.06)% to x= 1.37(.12)%. At least 
10 locations were measured on each crystal surface as a 
means of ensuring that the stoichiometry was homoge- 
neous. In all but one sample, the uniformity was within 
the statistical limitations of the measurement. For one 



sample it was necessary to examine 54 different locations 
covering the crystal surface to map out and take account 
of a small region, approximately 5% of the total surface, 
having an elevated Gd concentration of x= 1.50(.08)%. 
Measurements on the cut edges of a 100 [im thick sample 
confirmed that there was no gradient through the thick- 
ness of the sample. If the Gd concentration in the melt 
were to drift with time during the crystal growth, then a 
resulting Gd gradient in the crystal would be problematic 
for the following analysis. 

The broadband microwave absorption spectrum 
(0.5 GHz to 21 GHz) was measured using a novel spec- 
trometer based on a bolometric method of detection, de- 
scribed in detail elsewhereiSS In essence, the single crystal 
GdxYi_xBa2Cu306+y sample is exposed to a spatially 
uniform microwave magnetic field that is amplitude mod- 
ulated at low frequency. The corresponding temperature 
oscillations of the sample are detected synchronously, and 
absolute calibration of the thermal stage sensitivity al- 
lows us to infer the magnitude of the power absorbed 
by the sample. A critical ingredient of the method is 
a well-characterized normal-metal reference sample that 
is mounted on a second thermal stage and placed in a 
symmetric location in the microwave transmission line. 
This allows us to monitor the strength of the applied mi- 
crowave field which varies strongly with frequency due to 
standing waves in the microwave circuit. 

At microwave frequencies in the limit of local 
electrodynamics,^"'^ the experimentally measurable quan- 
tity is the surface impedance Zg which is related to the 
complex conductivity a ^ g\ — \ai via 

Z,EER,+iX, = .p^. (1) 

For a superconductor well-below measured at low fre- 
quency, the high superfluid density ensures that the re- 
sponse is mainly reactive (i.e. (72 S> cri) and Eq.^ sim- 
plifies to 

R,(c.,T) ~ \^lu?WT)Gr{uj,T), (2) 
X,(c.,T) ~ [^=^^coXiT). 

V 0-2 

(The appearance of A in both of the above expressions 
highlights its important role in determining the conduc- 
tivity from measurements of Rs and Xg.) The power 
absorption is determined by the surface resistance Rs = 
Re |Zs I according to 

Pabs = Rs J H^fdS, (3) 

where Hrf is the root-mean-square (rms) magnitude of 
the uniform, tangential magnetic field at the surface S of 
the sample. At low temperatures, absorption in a super- 
conductor is due to quasiparticles thermally excited from 
the condensate, and the study of the quasiparticle con- 
ductivity spectrum cri (w, T) has been the central focus 
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FIG. 2: Measurements of the temperature-dependent change 
in apparent magnetic penetration depth, AA(r), performed 
via cavity perturbation in a loop-gap resonator operating at 
954 MHz. A comparison between measurements on single 
crystal samples with and without 1% Gd substituted for Y 
demonstrates that only below a temperature of ~7 K does the 
contribution from the real part of the susceptibility x'i^jT) 
contribute to the measurement as an upturn in the apparent 
penetration depth. Above this temperature, the intrinsic lin- 
ear dependence of AA(T) of a clean d-wave superconductor is 
recovered. 



of other recent workiSiiS^ In the present case of magnetic 
impurities in a superconductor, the imaginary part of the 
magnetic susceptibility, x"j provides another mechanism 
for power absorption. Formally, the magnetic response 
can be included by a modification of the vacuum perme- 
ability fj,o in Eq. ^ to include a dimensionless complex 
susceptibility x, written as /io = /^oll + x' ^ ix")- 
For the case of a sufficiently low concentration of mag- 
netic impurities, x\x" ^ 1: ^^id the apparent surface 
impedance expressions can then be rewritten as 



(Ti(w,r) + 



x"(^,T) 



(4) 



Thus, a frequency scanned measurement of the power ab- 
sorption will provide an apparent surface resistance with 
two separate contributions: one from the quasiparticle 
conductivity (Ti(w,T) and one from the ESR spectrum 
contained in x"(w,T). 

Equation 0] reveals that a measurement of the surface 
reactance will contain a contribution from x'(cl;,T) in 
addition to probing A(r). This is evident in the com- 
parison of our measurements of AA(T) in pure and Gd 
doped YBa2Cu3 06.99, presented in Fig. 13 In fact, at the 
outset of this work we considered the possibility of fitting 
the curvature arising from x'{T) in a low frequency mea- 
surement of AA(T) as a means of extracting the absolute 
value of A. However, it was decided that this method 
was not practical as it lacked the necessary sensitivity 
and also required knowing x'(a;,T) a priori. A further 
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FIG. 3: Broadband microwave absorption measurement of a 
Gdo.oi Y0.99Ba2Cu3O6.9g single crystal sample at 1.3 K. The 
broad, smooth background is due to quasiparticle absorption 
and the sharp lines are due to ESR transitions. The smooth 
line is a polynomial fit to the background away from the ESR 
peaks, and is seen to be similar to a sample without Gd. The 
inset shows that the main ESR transition at 14.05 GHz is very 
well described by a Lorentzian lineshape, shown as a dashed 
line. 



complication to this approach is that it has been well- 
established that a variety of paramagnetic impurities can 
generate a low temperature upturn in a low frequency 
AA(r) measurement,^'^ which would produce systematic 
errors in the extraction of an absolute value of A. 

Fortunately, the broadband spectroscopy apparatus 
provides a means of accessing the ESR spectrum in 
much greater detail. Figure El shows a measurement of 
R°PP(u;) from 0.5 GHz to 21 GHz at 1.3 Kelvin for a 
Gd0.01Y0.99Ba2Cu3O6.99 sample. The quasiparticle ab- 
sorption spectrum is known to be a slowly varying, mono- 
tonically increasing function of frequency, and the sharp 
ESR absorption hues (Af < 0.5 GHz) are superposed. 
We have shown previously that the low temperature mi- 
crowave conductivity ai{u>,T) extracted from measure- 
ments of Rs{uj,T) can be interpreted as that of weak- 
limit quasiparticle scattering from impurities;^ and here 
we attribute the small changes in the spectrum of Fig. El 
to a slight increase in scattering due to the increased den- 
sity of out-of-plane defects in the Gd doped samples. In 
order to separate the ESR contribution from the quasi- 
particle conductivity, we fit to the smooth background 
away from the ESR peaks using a second order polyno- 
mial, shown as a solid line in Fig. El After subtract- 
ing the polynomial, the remaining quantity is equal to 
^fio^^Xx"{oj) (from Eq.0. In the subsequent section, we 
discuss the calculation of the quantity x"('^) from the 
effective spin Hamiltonian. With the theoretical curve 
for the susceptibility of one Gd spin in hand, we are able 
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to fit to the experimentally determined quantity with an 
overall multiplicative fit parameter A. 

YBa2Cu306+y single crystals grow naturally as 
platelets having broad a6-plane crystal faces and thin 
c-axis dimensions. All of our microwave measurements 
involve a uniform radio frequency magnetic field H^f ap- 
plied parallel to the broad face of the crystal in order 
to minimize the effects of demagnetization, as shown 
schematically in Fig. ^ It is also known that the crys- 
tal c-axis is the principal spin axis of the Gd ions, so 
that the ESR response is maximized in this configura- 
tion. The number of Gd spins that are exposed to the 
microwave field is governed by the effective penetration 
depth Ae// and, for our geometry, will always contain 
contributions from Aq or Ab and Ac. In the cuprates, 
the two-dimensionality of the Cu02 planes results in a 
large anisotropy in the penetration depth: Ac ^ Aa,Af, 
for all doping levels. In practice, this means that despite 
working with thin crystals having the in-plane dimension 
^a.b ^ lOfc, the c-axis contribution to X^ff can be large. 
In some cases it was also necessary to account for a small 
fraction t of the crystal which remained twinned (14% 
at most), introducing the third component of A into the 
measurement. Since A is small compared to any crystal 
dimension, it suffices to use a simple linearized relation 
for Ae//, 

^ lp{{l-t)\a+t\b)+lc\c 

Mf = 7 (5) 

where £p is the in-plane crystal dimension parallel to the 
flow of diamagnetic screening currents (the geometry is 
that of Fig.^. Each component of A can be isolated by 
combining measurements performed on the same crys- 
tal but with its geometry altered in a controlled manner 
as attained by rotating and cleaving the crystal. This 
matter is discussed in detail in Appendix IXI 



III. THE CRYSTAL FIELD HAMILTONIAN 

In this section, the way in which the magnetic penetra- 
tion depth is obtained from our measurements of the ESR 
absorption spectrum is explained. This amounts to the 
treatment of a dilute random array of Gd'^"'" ions, each 
having an electron spin oi S = 7/2, that contribute to 
the measured quantity ^fi^Xx" {i^)- We point out that in 
this work the small splitting due to the Gd isotopes with 
non-zero nuclear spin can be ignored.-^ The spectrum 
is a result of the energy level configuration determined 
by the splitting of the degenerate Gd spin levels by the 
crystalline field. Because the spin system is dilute, it 
can be described by a non-interacting single-spin effec- 
tive crystal field (CF) Hamiltonian. The application of a 
microwave field, oriented perpendicular to the principal 
spin axis of the system, is a time dependent perturbation 
that induces transitions between the spin levels. When 
the frequency of the applied field is tuned to the splitting 



between two levels, a maximum in the ESR absorption is 
observed. The intensity of the transition is found through 
the application of Fermi's golden rule for the two levels 
that define the transition, rrii and m-,-, which allows the 
transition's contribution to the susceptibility to be writ- 
ten as 

^IM = A^oJlAi^.f '^(c^ - [E, - E.,]/h)[N, - N,]. (6) 

In this expression, the eigenstates rrii and rrij correspond 
to the energies Ei and Ej . The level occupation number 
is Ni = iVo-Z"^ exp(— /^E'i), where A^o is the number of 
spins per unit volume, Z is the partition function and 
[3 = l/ksT is the inverse temperature. The time depen- 
dent perturbation induced by the microwave field, Hrf , is 
proportional to the Sx or Sy operator and the matrix el- 
ement between the two states is /i^ — g^B{'mj\Sx/y\mi) , 
where g is the Lande factor and /is is the Bohr magneton. 

In real materials the ESR transitions are broadened 
by various relaxation processes such as spin-spin interac- 
tions, lattice disorder and spin-lattice interactions. This 
broadening, however, does not change the overall inten- 
sity of the ESR line, but merely replaces the delta func- 
tion by some line shape function (in our case it is close 
to Lorentzian) whose integrated intensity is unity. In the 
initial stages of this work we examined optimally doped 
(y=0.93) crystals having three different concentrations 
of Gd ions; x~ 0.5%, x~ 1%, and x~ 3%. From these 
measurements, we were able to conclude that the width 
of the Lorentzian-shaped ESR peaks {5i = 0.35 GHz for 
x= 1%) scaled approximately linearly with the nominal 
Gd concentration x, but with a substantial x=0 inter- 
cept of about 0.30 GHz. Although we did not perform 
the detailed EPMA analysis to measure x of the nomi- 
nally 0.5% and 3% samples, it seems clear that at x= 1% 
there is a dominant contribution to the broadening that 
is concentration independent. Since a dilute spin-spin 
interaction is expected to provide a relaxation rate that 
is linearly proportional to the concentration's our re- 
sults suggest that some other process is responsible for 
much of the line broadening. The line width was also 
measured in the optimally doped samples at two differ- 
ent temperatures as shown in Fig. ^ The line widths 
decrease slightly upon increasing the temperature from 
1 K to 3 K, indicating that spin-lattice relaxation is not 
likely to be the dominant process since typical relaxation 
mechanisms would be strongly temperature dependent. 
It remains unclear what the main spin relaxation mech- 
anism responsible for the ESR line width is here. A final 
point is that no power dependence of the spectra were 
ever observed verifying that saturation effects were neg- 
ligible and that the assumption of thermal equilibrium 
implicit in the Boltzmann factors used in Eq. is well- 
founded. 

For the purpose of extracting A, the quantity of in- 
terest is the effective number of spins exposed to the 
field Hrf. This requires a knowledge of the energy lev- 
els, their population and the corresponding matrix ele- 
ments. This in turn requires knowing the Hamiltonian 



6 



o 3% Gd, 1 .3 Kelvin 
A 3% Gd, 3.0 Kelvin 



c 

3 




11 



12 



13 14 
Frequency (GHz) 



15 



16 



FIG. 4: The temperature dependence of the ESR hne 
width for 1.3 K and 3.0 K of an optimally doped 
GdxYi_xBa2Cu306+y sample with x=.03, normalized by the 
height of the 1.3 K peak. The G-band has line widths of 0.595 
and 0.571 GHz and the F-band has line widths of 0.649 and 
0.636 GHz at 1.3 K and 3 K, respectively. The line widths in 
our measurements of 1% Gd samples presented in this work 
are widths of 0.25-0.46 GHz at 1.3 K. 



±1/2 
±3/2 



I ±5/2 

LU 



±7/2 











> 









rO rO 

£»2'^4 



» 

02,04, lSf^,tS 



,0 
>6' 



FIG. 5: The energy levels resulting from the effective crystal 
field Hamiltonian. The levels, from left to right, correspond 
to: a diagonal Hamiltonian, H oc O2 = 3S^ — S{S + 1), a 
tetragonal Hamiltonian with some off-diagonal level mixing 
and finally an orthorhombic Hamiltonian with O2 = (<S'+ -I- 
S'L)/2. The CF parameters shown here are typical for the 
samples we have measured, except for B2 which has been 
exaggerated in order to make its effect more visible. 



of the system. The Gd atom has the electron configu- 
ration of [Xe]4/^5c?^6s^ and ionization number of three 
in GdxYi_xBa2Cu306+y. This implies that the outer 
shell of the ion remains exactly half filled (seven / elec- 
trons in fourteen states) and the system behaves like a 
7/2 spin with no mixing with other S + L multiplets 
due to Hund's rule. The construction of an effective 
Hamiltonian for a spin system is based on the expan- 
sion of the crystal field around the magnetic ion in terms 
of Stevens operators. These are the operator analogues 
of the spherical harmonics and are functions of S, Sz 
and S±, given, for example, by Abragam and Bleaneyj^^ 
Due to the lattice symmetries and the finite spin of the 
magnetic ion, the expansion reduces to a finite sum of 
operators H = ^ B^O^ where the are numerical 
coefficients called the crystal field parameters. The high- 
est order of Stevens operator allowed for a spin of 
7/2 in a crystal field is 6 since the matrix elements of 
the Stevens operators between two spin states are zero 
unless p < 2S. The order is restricted to even values by 
the required time reversal symmetry of the Hamiltonian. 
Furthermore, the a6-plane anisotropy is so small that for 
most purposes one can assume 7r/2 rotation symmetry 
about the c-axis; orthorhombicity is most significant in 
the overdoped samples. In the case of tetragonal sym- 
metry, the relevant operators in order of significance are 
0^,0l,0j, and O^. The sixth order terms are dif- 
ficult to resolve in our measured spectra and we simply 
adopt those of previous authorsi2Si2L2& When the tetrag- 



onal ordering is slightly distorted, as will be discussed in 
the next section, three more Stevens operators can con- 
tribute to the Hamiltonian. For the present work, the 
largest and only significant one of them is 0| . 

The spin 7/2 system discussed herein has four dou- 
bly degenerate energy levels with four eigenstates that 
are roughly the ±Sz eigenstates of the Sz operator. The 
three allowed ESR transitions correspond to Sz •S'zibl. 
Small off-diagonal terms in the Hamiltonian, such as the 
and , induce small level mixing that allow the oth- 
erwise forbidden transitions corresponding to ±5/2 — s- 
±1/2, ±7/2 ^ ±1/2 and ±7/2 ^ ±3/2 to occur, al- 
beit with relatively low intensities. The resulting energy 
levels for different CF Hamiltonians are shown in Fig. |S1 
Due to the oxygen ordering in chains, which we discuss 
in detail in the next section, most samples have more 
than one inequivalent site for the Gd ions. Therefore, 
two or three sets of CF parameters are needed to fit the 
spectrum. Following Janossy et aLrSi2L2& we label the 
spectra corresponding to a particular site by a band index 
using the notation shown in Fig. |S1 

We begin our fitting procedure by first constructing 
an effective spin Hamiltonian using the measured CF pa- 
rameters reported by Janossy et alm^ and then adjusting 
the parameters to best fit our data. We focus on the 
main transition (±7/2 ±5/2 for each band) because 
it is the strongest and consequently least susceptible to 
experimental uncertainty. As a final step, a multiplica- 
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tive factor representing the effective number of spins that 
have participated in the ESR process is then used to scale 
the overall amplitude of the model spectrum to match the 
measured data. This number is the effective penetration 
depth Xeff, combined with the measured Gd concentra- 
tion X in the crystal. 



IV. DEPENDENCE OF ESR SPECTRA ON 
OXYGEN CONFIGURATION 

The presence of a chain layer breaks the tetragonal 
symmetry and the Stevens operator expansion of the 
foregoing section is then no longer exact. However, the 
chains are far from the Gd sites and this effect can 
be treated as a perturbation, leading to line broaden- 
ing or line splitting. Here we will consider line split- 
ting associated with distinct Gd environments and sim- 
ply handle the line width as a parameter in a Lorentzian 
fit. Different types of Gd sites are encountered because 
hole doping in this system is controlled by manipulating 
the oxygen content and ordering of the CuOy chains in 
GdxYi_xBa2Cu306+y The oxygen content is set by an- 
nealing in controlled oxygen partial pressure at high tem- 
peratures and then the oxygen ions tend to organize into 
lengths of CuO chain fragments which are able to pro- 
mote holes to the Cu02 planesi^ The strong tendency 
to form chain fragments means that there are 7 prob- 
able Gd environments corresponding to anything from 
to 4 nearest neighbor chains. These different crystal- 
lographic environments have been identified in ESR ex- 
periments on magnetically aligned GdxYi„xBa2Cu306+y 
powdersSSiSLSS and are illustrated in Fig. 

The chains also tend to form ordered periodic superlat- 
tices consisting of arrangements of full and empty CuOy 
chains. The structures can be particularly well-ordered 
for special values of y. In the case of full oxygen doping, 
the CuO chains are nearly completely full with every Gd 
ion having four nearest neighbor chains. This results in 
the simplest of the spectra observed in the five crystals 
studied. It involves one primary set of CF parameters, 
denoted the G-band by Janossey et ai, and generates 
three dominant ESR transitions. A secondary contribu- 
tion of much lower intensity is also present in these mea- 
surements, resulting from a configuration where one of 
the four chains is missing oxygen, denoted the F-band. In 
the overdoped sample, the largest transition (7/2 5/2) 
of this latter band is barely discernable, but as the oxygen 
concentration is reduced to optimal doping, the intensity 
of the F-band builds, as seen in Fig. [7| 

At lower doping, the next best-ordered phase occurs 
near y « 0.5, where the chains form an ortho-II struc- 
ture of alternating full and empty chains.^® Other ordered 
phases are: ortho-III with FEE (full, full, empty) chains 
and ortho-VIII with FFEFFEFE chains. In each of these 
phases there is more than one possible chain configura- 
tion around the Gd ions. Having detailed spectra at each 
doping allows us to identify the different bands and fit 



each spectrum with a set of crystal field parameters for 
each band. The fitted x" spectra are shown in Fig. [7| 
In each spectrum, the relative intensities of the different 
bands is a measure of the relative number of Gd ions in 
each configuration. This information can be used as a 
measure of the amount of oxygen in the sample: 

n n 

where n is the number of full nearest neighbor CuO 
chains in the band and /„ is the observed relative inten- 
sity of the line associated with a particular band. The 
oxygen content is given by 6 -I- y. The ESR measured 
oxygen content versus the chemically measured content 
is presented in Fig. |S1 Note that the ESR points consis- 
tently underestimate the amount of oxygen in the crys- 
tals. This is likely due to the inefficiency of isolated oxy- 
gens in empty chains in promoting a hole into the Cu02 
planes. The above oxygen content analysis assumes max- 
imal length of chains, i.e. all chains are either full or 
missing with no isolated oxygens or vacancies. Since the 
Gd ions are likely to be sensing the presence of chains 
through their influence on the charge distribution on the 
much nearer Cu02 planes, an isolated oxygen ion hardly 
affects the crystal field environment since it does not pro- 
mote a hole into the planes. Thus a Gd ion near an iso- 
lated oxygen in an otherwise empty chain will experience 
the crystal field of an empty chain. This will lead to an 
underestimate of the oxygen content, but will not affect 
the overall counting of Gd ions. 

It is also interesting to note that the ESR analysis pro- 
vides a good measure of the chain disorder. For exam- 
ple, perfect three dimensional ortho-II ordering should 
produce only one band, corresponding to two full chains, 
the D-band. However, our ortho-II spectrum displays 
two additional bands (the B- and F- bands) resulting 
from imperfections in the ortho-II ordering^ These de- 
viations can be interpreted as ortho-II phase boundaries. 
Our ortho-III and ortho-VIII spectra are consistent with 
long range ortho-III and ortho-VIII order along the ah 
plane but random stacking along the c-direction. 



V. RESULTS 

We have performed detailed measurements of the ZF 
ESR absorption spectrum in high quality samples of 
Gdo.oiYo.g9Ba2Cu306+y at five different values of y. As 
noted above, four of these are chosen to be in ranges with 
particularly well-ordered CuO chain structures, denoted 
by the periodicity of their chain superlattices (ortho- 
I, II, HI, VIII). The optimally doped sample, with the max- 
imum attainable Tc of 93 K (y=0.93) is also in the ortho-I 
phase, but has a substantial number of oxygen vacan- 
cies on the chains. The ESR spectrum at each doping 
was fitted with a spectrum generated by an effective CF 
Hamiltonian, with initial guesses for the CF parameters 
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FIG. 6: Identification of the ESR bands (A,B,D,F,G) attributed to a variety of oxygen configurations, taken from S. Pekker 
et al^ Fi gures represent cross-sections through the center of a YBCO unit ceU (with Gd located on the Y site) paraUel to the 
ac-plane. Open circles represent doped holes in the Cu02 planes, black circles represent oxygen in the CuO chains and diamonds 
represent oxygen vacancies in the CuO chains. The central circle is the Gd ion. Note that the three D-band configurations all 
result in a single ESR band. 
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FIG. 7: Measured ESR spectra for five dopings with screening currents running along the a direction, shown together with 
the spectrum calculated using the CF Hamiltonian. The calculated curve is scaled by an overall multiplicative factor of Ae// as 
explained in the text. The various spectra were offset by 0.1 from each other and the ortho-I spectrum was scaled down by a 
factor of four. 



taken from Rockenbauer et alm^ Subsequently, the CF 
parameters were adjusted to best describe our measured 
spectra, and the best fit values are presented in Table IVl 
The evolution of the dominant CF parameter (S^) as a 
function of the ESR measured oxygen content is shown in 
Fig. 1^ The values and systematic trends of our fit values 
compare well to those of the conventional ESR studies 
on powders i^^i^LSS, 

The ESR spectrum is a sensitive probe of the crystal- 



lographic structure and the crystal symmetry in partic- 
ular. This high sensitivity is demonstrated at its best 
in the spectrum of our overdoped sample. This is the 
most ordered and most orthorhombic phase, with four 
full CuO chains around virtually every Gd ion. The lines 
are very sharp and the main band, G, corresponds to 
96% of the total spectral weight. As seen in Fig.^] the 
orthorhombicity of the crystal is manifested as an addi- 
tional line in the spectrum that appears only when the 
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FIG. 8: Ordered oxygen content as measured by ESR vs. the 
total oxygen content of the sample. The solid circles are taken 
from Pekker et al?^ (ytotai determined by weighing the sam- 
ple) and the open diamonds are the values obtained from the 
present work (ytotai determined by the annealing temperature 
and oxygen partial pressure). 



measurement is performed with screening currents flow- 
ing along the 6-direction. The spectra for both directions 
can be fit with a single set of CF parameters by introduc- 
ing an orthorhombic term in the CF Hamiltonian. This 
term is i?|0|, where _B| is the measured coefficient and 
Ol = {Si + 81)12 = {Si - Sl)/2. The inclusion of 
this term allows the otherwise forbidden transitions of 
±5/2 ±1/2 only when the magnetic field is applied 
in the a-direction, and reveals the slight orthorhombic- 
ity of the crystal (the difference between the d and h- 
dimensions of the unit cell is less than 2%). 

The measured ortho-Ill spectrum is indicative of some 
disorder in the chain structure. First, the substantial 
amount of spectral weight in the F-band, which corre- 
sponds to three full chains, indicates that even if the 
ordering is perfect in the planes, the stacking along the 
c-direction is disordered, i.e. that the ortho-III struc- 
ture may shift between adjacent chain layers. Second, 
we find that each ESR peak is composed of two closely 
spaced lines that together produce a non-Lorentzian 
shape. Since the measured ESR oxygen level is consis- 
tent with the other crystals (see ytotai = 0.77 in Fig.|SJ), 
we are confident in the band identification and do not 
interpret the minor splitting as a change of chain con- 
figuration. The splitting implies a small variation of 
CF parameters within the same band. This might be 
the result of inhomogeneity in the oxygen concentration, 
which is consistent with a measured transition temper- 
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FIG. 9: The oxygen doping evolution of the main CF param- 
eter, B2 in MHz for the different bands. The open symbols 
are taken from Rockenbauer et al^ and the stars are the ZF 
ESR best fit values. For the ortho-III sample, where more 
than one set of CF parameters was found for each line, we 
show the values closest to previously the measured ones. 



ature broadening ATc ~1 K. However, we do not think 
that the sample is macroscopically phase separated since 
the line shapes were not affected by cleaving the sample 
into smaller pieces. 

Our ortho-II spectra suggest fairly good ordering 
within the chain layers with a somewhat shorter correla- 
tion length in the c-direction. This is in agreement with 
X-ray analysisii^ The extraction of Ac in this sample was 
done by comparing two measurements in the a-direction, 
where the first measurement was done on the whole sam- 
ple and the second one on two pieces resulting from a 
cleave along the 6-direction. The effective penetration 
depth, Ae// was observed to increase by about 210 nm due 
to the introduction of two more sample 6c-faces, where 
currents flow along the c-direction. 

We have obtained enough data in order to reliably ex- 
tract the absolute values of the penetration depth in all 
three crystallographic directions for three of the five dop- 
ing levels studied: fully doped, ortho-II and ortho-III. 
Our measurements at the remaining two doping levels 
were made using samples that were not amenable to the 
cleaving procedure required to extract absolute values of 
A. Our final results for A were derived as described in 
the previous sections with the emphasis on performing 
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FIG. 10: The main ESR line ±5/2 ±7/2 in ortho-I, ortho-III and ortho-II. In each figure we overlay three or four 
measurements of the same crystal in the a and 6-directions, before and after cleaving the sample. Each graph has been scaled 
by the effective penetration depth X^ff of the particular measurement so that only x" is presented. The changes in X^ff are due 
to the different aspect ratios of the crystal before and after the cleave and the different orientations of measurements. 



enough measurements on the same sample in order to 
eliminate sample-dependent effects and to overdetermine 
the values of A for error control. The penetration depth 
results are presented in Table |nj The error estimates 
account for uncertainties arising from measurements of 
the sample dimensions, calibration of the microwave ab- 
sorption experiment?^, and the estimation of the ESR 
fit parameters. Since the extraction of the penetration 
depth in the three crystallographic directions depends 
sensitively on the sample's dimensions, it is important 
to work with crystals that have smooth parallel well de- 
fined faces whose area is easy to measure. This was most 
easily achieved by choosing a single sample having a nice 
platelet geometry for each doping. If we restrict our anal- 
ysis to measurements made on the same crystal where 
the geometry of the sample can be altered in a controlled 
fashion by cleaving, we obtain a high degree of internal 
consistency between different measurements. In all cases 
we performed an extra measurement that overdetermines 
the penetration depth values and found agreement rang- 
ing from 0.5% to 6%. In the case of the ortho-II doping 
we measured a second sample having a 30% higher Gd 
concentration, and the results agreed to within 12% of 
the first sample. The high reproducibility of our mea- 
surements is demonstrated in Fig. 1101 where we present 
data for the same crystal before and after cleaving the 
sample. Despite changing the aspect ratio of the crystal 
by approximately 50%, the spectra are identical, modulo 
the overall multiplicative factor of Ae// . 

VI. DISCUSSION 

The results presented here have impact in two areas. 
First, we have provided detailed information about the 
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FIG. 11: The oxygen-overdoped Gd.01Y.99Ba2Cu3O6.993 
sample ESR absorption spectra in the a and 6-directions. 
The 6-axis data were offset by 25 nm for clarity. The or- 
thorhombicity of the crystal lattice is apparent in the appear- 
ance of an additional ESR line at 6.9 GHz when the screen- 
ing currents flow in the 6-direction. This line corresponds 
to the ±1/2 ^ ±5/2 transition that is allowed along the b- 
direction, generated by the orthorhombic Stevens operator 
B2 = {S^ — Sy)/2. This same line appears with suppressed 
intensity in the a-direction due to a small fraction of twinning 
in the sample. The observed spectrum indicates a predom- 
inance of the G-band oxygen configuration, as expected for 
a fully doped sample, with a 4% contribution from the F- 
band. For comparison, note that the F-band contribution in 
the optimally doped sample shown in Fig.|7]is much larger as 
a result of the lower oxygen content. 



11 



Sample 
(Gd%) 


Band 


Intensity 


(MHz) 


Bl 
(MHz) 


Bl 
(MHz) 


Bt 
(MHz) 


7/2 ^ 5/2 
(GHz) 


3/2 ^ 5/2 
(GHz) 


1/2 ^ 3/2 
(GHz) 


ortho-I 
(1.08 ±0.06) 


G-band 


96% 


-587 (-598) 


26.7 (-) 


-3.1 (-3.1) 


12.5 (13.2) 


14.0 


5.8 


1.12 


F-band 


4% 


-526 (-541) 


26.7 (-) 


-3.1 (-3.1) 


12.5 (17.2) 


12.9 


5.0 


0.9 


optimal 
(1.28 ±0.02) 


G-band 


75% 


-590(-598) 





-3.1 (-3.1) 


12.7 (13.2) 


14.04 


5.86 


1.05 


F-band 


25% 


-533(-541) 





-3.1 (-3.1) 


12.5 (17.2) 


13.02 


5.2 




ortho-III 
(1.02 ±0.05) 


G-band 


21% 


-617,-604,-585(-598) 


30 (-) 


-3.0 (-3.1) 


12.5 (13.1) 


14.4,14.1 


6.13,6.03 


1.4,1.3 


F-band 


56% 


-533,-517 (-541) 


6.7 (-) 


-3.5 (-3.1) 


17.1 (17.2) 


13.3,13.0 


5.4,5.45 


0.2,0.3 


D-band 


23% 


-481,-475 (-497) 


26.7 (-) 


-3.1 (-3.1) 


14.6 (14.7) 


12.1,11.9 


4.7,4.6 


0.6 


ortho-VIII 
(1.38 ±0.1) 


G-band 


2% 


-630 (-598) 


1.6 (-) 


-2.5 (-3.1) 


9 (13.1) 


14.2 


6.3 


1.9 


F-band 


36% 


-549 (-541) 


1.6 (-) 


-3.0 (-3.1) 


12.5 (17.2) 


13.2 


5.5 


0.8 


D-band 


62% 


-497 (-497) 


1.6 (-) 


-3.1 (-3.1) 


12.5 (14.7) 


12.3 


4.9 


0.8 


ortho-II 
(1.08 ±0.06) 


F-band 


4% 


-560 (-541) 


5 (-) 


-3.1 (-3.1) 


17.2 (17.2) 


13.4 


5.8 


0.8 


D-band 


81% 


-502 (-497) 


5 (-) 


-3.2 (-3.1) 


14.7 (14.7) 


12.5 


5.0 


0.4 


B-band 


15% 


-437 (-445) 


5 (-) 


-3.2 (-3.1) 


15.5 (17.2) 


11.2 


4.5 


0.2 



TABLE I: Crystal field parameters for the various ESR bands for different oxygen contents of Gd.oi Y,99Ba2Cu306+y, given in 
MHz. The Gd% in brackets indicates the relative substitution of Gd for Y. The CF parameters measured by Rockenbauer et 
al.^^ in conventional ESR on powder samples are given in brackets, and at each doping we have taken values from the nearest 
doping measured by these authors. The values of ,65 and Bg could not be resolved from our spectra and we simply adopt the 
values of Rockenbauer et al. The transition frequencies, shown in GHz, were found either directly from the observed lines or by 
the fit Hamiltonian (when the lines were unobservable) . The ortho-HI spectrum required more than one set of CF parameters 
for each band due to further line splitting of unknown origin, which may be related to oxygen inhomogeneity. 



Crystal Structure 


rc(K) 


Oxygen Content 


Aa(nm) 


A6(nm) 




Ac(nm) 


Ortho-I 


89 


6.995 


103 ±8 


80 ±5 


1.29 ±0.07 


635 ± 50 


Ortho-III 


75 


6.77 


135 ± 13 


116 ± 12 


1.16 ±0.12 


2068 ± 200 


Ortho-II 


56 


6.52 


202 ± 22 


140 ± 28 


1.44 ±0.26 


7500 ± 480 



TABLE II: Experimental values of the anisotropic magnetic penetration depth extracted from four or more different mea- 
surements of \eff on the same crystal, with the microwave field applied in different directions and for different sample aspect 
ratios achieved by cleaving. The uncertainty in the ansisotropy is reduced since some of the systematic contributions to the 
uncertainty in absolute value cancel when comparing measurements on the same sample. 



GdxYi_xBa2Cu306-i-y system, namely the crystal field 
parameters at the yttrium site and the way in which 
oxygen orders into CuO chains. We have shown that 
our crystal field parameters agree very well with those 
obtained by previous authors using the completely differ- 
ent method of high field ESR on powdered samples. Our 
results also support the conclusions of the same work 
which found that isolated oxygens in vacant chains do 
not influence the CF conflguration around the magnetic 
ion and therefore cannot be distinguished from a vacant 
chain via ESR spectroscopy. This is in accord with the 
fact that isolated oxygens do not promote holes from the 
copper-oxygen planesi^S 

The principal result of this work is a set of new 
values for the penetration depth in all three crystal- 
lographic directions for three oxygen-ordered phases of 
YBa2Cu306-i-y which are summarized in Table II. The c- 
axis penetration depth is large and increases very rapidly 
with decreasing doping, as observed previously in in- 
frared measurements of Ac by Homes et ali^ The in-plane 
measurements are in accord with the results of Sonier et 



4,30 -^j^o performed fiSR measurements on mosaics of 
high-purity, dctwinncd single crystals at two other oxy- 
gen dopings. The muon measurements are unable to di- 
rectly determine anisotropics, but they compare well with 
the geometric mean of the a5-plane values reported here. 
A sample of YBa2Cu3O6.60 with Tc = 59 K was found to 
have Xab = 170 nm^S very close to the geometric mean 
VAaAfo = 168±26 nm of the Tc — 56 K sample studied 
here. An optimally doped sample of YBa2Cu3 06.95 with 
Tc = 93 K had Xab = 112 nm,^ which lies between the 
means of the Tc — 75 K sample (125±12 nm) and the 
overdoped Tc = 89 K sample (91±6 nm) studied here. 
The agreement reflects the particular care taken in these 
fiSR measurements to cover a wide range of applied mag- 
netic fields in order to ascertain the low field limiting 
values and thus minimize the non-linear, non-local, and 
other effects that can arise in the vortex state. It must 
also be noted that the fj,SR measurements hinge on a 
detailed model of the vortex lattice that provides a fit to 
the field distribution detected by the muons. 

These absolute values of penetration depth obtained 
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in single crystals are smaller than many measurements 
found in the existing literature on YBa2Cu306+y. Mag- 
netization studies of aligned powders by Panagopoulos 
et ai^ gave 140 nm for Tc^92, 210 nm for Tc=66 and 
280 nm for Tc—56 (all ±25%), which arc longer than the 
values presented here. The origin of this discrepancy is 
not clear, but might be due to assumptions made in the 
analysis of the powder data or problems with the sur- 
faces of grains embedded in epoxy. Our absolute values 
of Aq and A;, are also smaller than those obtained by far 
infrared measurements near optimal dopingii (Aa=160 
nm, A{,=100nm). However, the infrared measurements 
did point out the importance of the in-plane anisotropy in 
these materials, an anisotropy that becomes very large in 
the YBa2Cu3 6.99 sample studied here. This anisotropy 
has been attributed to the presence of a nearly one- 
dimensional Fermi sheet derived mainly from the bands 
associated with the CuO chainsi^ Although jiSR mea- 
surements cannot directly measure this anisotropy, Tal- 
lon et alm^ also inferred such a contribution by noting 
the very large increase in muon depolarization rates as 
the chain oxygen sites become filled near YBa2Cu307. 
The interpretation of the polycrystalline data suggested 
values of Xa = 155 nm and Xb = 80 nm for samples 
with fully doped CuO chains, which overestimates both 
the overall magnitude and the anisotropy of the in-planc 
penetration depth. However, the basic picture of chain- 
driven anisotropy is supported by our new microwave 
measurements, which also show substantial anisotropy in 
the ortho-II ordered sample, but rather less in the ortho- 
ni sample, a sensible trend since the ortho-III ordering 
is much poorer, leading to more fragmented chains. 

Much of the work on the doping-dependence of Xab has 
relied on muon spin relaxation measurements on poly- 
crystalline samples. In these measurements, the field in- 
homogeneity in the vortex state gives rise to dephasing 
of the precessing spins of implanted muons. A Gaussian 
fit is often used to extract a relaxation rate a which is 
deemed proportional to 1/A^. As long as the material is 
nearly two dimensional, the constant of proportionality 
can be calculated to be u = 7.09 x 10'*A~^, with a in /is^^ 
and A in nmi^^iM In the YBa2Cu306+y system, this very 
simple treatment of muon data gives values of Aaf, that 
are typically 20% or more larger than the single crystal 
data reported herei ^i^'^i'^^i^^ The resulting underestimate 
of the superfluid density, together with the problematic 
a6-plane anisotropy, mean that the single crystal pene- 
tration depths reported here can offer a clearer picture 
of the doping dependence of the superfluid density in the 
Cu02 planes. 

FigurelT^displavs Tc versus l/X^, a plot first suggested 
by Uemura et alM A long-accepted result has been a lin- 
ear relationship between these two quantities at low dop- 
ing, followed by a plateau at higher doping that varies 
from one cuprate system to another, inferred mainly from 
liSR measurements on ceramicsi^ The data in Fig. IT^ 
confirms the overall feature that 1/A^ increases with Tc- 
However, the linear regime at low doping, which implied a 
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FIG. 12: Measurements of the in-plane anisotropic absolute 
penetration depth A in YBa2Cu306+y for y=0.52, y=0.77 and 
y=0.995 plotted as Tc versus superfluid phase stifi'ness (oc 
l/A'^). Also shown are the fj,SR measurements of Sonier et al. 
on single crystal mosaics. The data do not support a case for 
scaling of Tc oc 1/A^ in this doping range above Tc = 56 K. 



superfluid phase stiffness proportional to Tc on the under- 
doped side of the cuprate phase diagram, is not supported 
by the new data presented here, at least for Tc > 56 K. 
The measurements on YBa2Cu3 06.993, which is slightly 
past optimal doping, are certainly up in the "plateau" 
regime for this material and would not be expected to 
follow Tc oc 1/A^. However, even for the lower two dop- 
ings with Tc = 56 K and Tc = 75 K the a-axis phase 
stiffness, which avoids contributions associated with the 
CuO chains, is still far from a relationship of the form 
Tc oc 1/A^. The plot instead suggests a sublinear depen- 
dence of Tc on 1/Aq, since Tc must fall to zero when 1/A^ 
does. Further data at low doping will be needed to clarify 
the relationship over the entire doping range. Interest- 
ingly, the latest muon spin relaxation studies by Tallon et 
al^ also suggest a sublinear relationship, although the 
values of A in that study are significantly larger than the 
single crystal data presented here. In the past, Tc oc 1/A^ 
has been seen as evidence that phase fluctuations play 
a central role in setting Tc on the underdoped side of 
the cuprate phase diagram. Indeed, the relatively low 
values of phase stiffness in the cuprates, plus the ob- 
servation of critical fluctuations near the superconduct- 
ing transistiou j^i^ do indicate that phase fluctuations 
play a role in determining Tc- However, the deviation 
from a linear relationship seen in Fig. E| suggests that 
other factors must also contribute to the critical temper- 
ature. The obvious candidate is thermal excitation of 
nodal quasiparticles^Si^^*^ which give rise to the linear 
temperature dependence of A(r) and rapidly deplete the 
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superfluid density as temperature rises. However, the 
central puzzle regarding the correlation between super- 
fluid density and Tc still remains in such a scenario. As 
the doping decreases, the zero temperature value of the 
superfluid density becomes much smaller, so that it be- 
comes easy for quasiparticle excitations, in addition to 
fluctuations, to drive the material normal at a lower Tc. 
The mystery still lies in understanding why the superfluid 
density becomes so small with decreasing hole doping. 
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APPENDIX A: EXTRACTING ANISOTROPIC A 
VALUES FROM MULTIPLE MEASUREMENTS 

YBa2Cu306+y single crystals have strong natural 
cleave planes perpendicular to both the [100] and [010] 
directions. Using this fact, previous studies of the tem- 
perature dependence of A have been successful at separat- 
ing the out-of-plane response from the in-plane response 
by measuring a sample, cleaving it N times, and then re- 
peating the measurement with all of the pieces together. 
The second measurement contains the same in-plane con- 



tribution, but has the c-axis contribution enhanced by a 
factor of iV.^^ For practical reasons, in the present work 
it was not possible to measure all of the crystal frag- 
ments together following cleaving in all cases, instead we 
successively measured each crystal after reducing the in- 
plane area. The c-axis thickness of the crystal was un- 
changed for all measurements, but the measurement of tc 
retains the largest relative uncertainty because £c ^ tab- 
Since Ac is also rather large, care must be taken when 
combining measurements in order to minimize the im- 
pact of the large relative uncertainty in the quantity ic^c 
in Eq. |21 The method we use here is to combine mea- 
surements of the same crystal in both a-axis and 6-axis 
directions before and after cleaving where the crystal's 
aspect ratio has changed by approximately 50%, provid- 
ing an overdetermination of the three unknown A values. 
For example, before and after cleaving the sample we 
measure 

^ a{{l - t)Xa + t\b) + cXc 
a + c 

_ bitXg + (1 - t)\b) + cAc 
which sum to give 

A second measurement with new dimensions t'^ and fj^ 
provide a second expression having the same form as 
Eg. lAll and subtracting the two eliminates Aq and Ah, 
leaving only Ac in terms of measured quantities. The 
values of Aa and Ab are then calculated from the above 
expressions. 
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